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Abstract 

We introduce an asymmetric classical Ginzburg-Landau model in a bounded inter- 
val, and study its dynamical behavior when perturbed by weak spatiotemporal noise. 
The Kramers escape rate from a locally stable state is computed as a function of the 
interval length. An asymptotically sharp second-order phase transition in activation 
behavior, with corresponding critical behavior of the rate prefactor, occurs at a crit- 
ical length similar to what is observed in symmetric models. The weak-noise exit 
time asymptotics, to both leading and subdominant orders, are analyzed at all interval 
lengthscales. The divergence of the prefactor as the critical length is approached is 
discussed in terms of a crossover from non-Arrhenius to Arrhenius behavior as noise 
intensity decreases. More general models without symmetry are observed to display 
similar behavior, suggesting that the presence of a "phase transition" in escape behav- 
ior is a robust and widespread phenomenon. 

KEY WORDS: Fokker-Planck equation, non-Arrhenius behavior, stochastic escape 
problem, stochastic exit problem, stochastically perturbed dynamical systems, spatiotem- 
poral noise, droplet nucleation, fluctuation determinant, instanton, false vacuum, stochastic 
Ginzburg-Landau models. 

1 Introduction 

Noise-induced transitions between locally stable states of spatially extended systems are re- 
sponsible for a wide range of physical phenomena^-'^l. In classical systems, where the noise 
*Partially supported by the National Science Foundation under grant PHY-0099484. 



1 



is typically (but not necessarily) of thermal origin, such phenomena include homogeneous 
nucleation of one phase inside another'^l , micromagnetic domain reversal'"^' ^' ^1 , pattern nu- 
cleation in electroconvection'^l and other non-equilibrium systems'^', transitions in hydrogen- 
bonded ferroelectrics'^1 , dislocation motion across Peierls barriers'^', instabilities of metaUic 
nanowires'-*^^!, and others. In quantum systems, the problem of tunneling between metastable 
states is formally similar, and problems of interest include decay of the false vacuumt-*^-*^! and 
metastable states in generall-*^^! , anomalous particle productiont-*^"^!, and others. 

The modern approach to these problems, beginning with the work of Langer on classical 
systems!^! and Coleman and Callan on quantum systemsl-*^-*^!, considered systems of infinite 
spatial extent (for a review, see Schulmant-*^^!). In certain systems, however, finite size may 
lead to important modifications, and in some instances qualitatively new behavior. Ap- 
proaches to noise-induced transitions between stable states in finite systems modelled by 
nonlinear field equations have been investigated by a number of authors'-*^^' 

In a recent paper'^^l, Maier and Stein studied the effects of weak white noise on a 
bistable classical system of finite size whose zero-noise dynamics are governed by a symmet- 
ric Ginzburg-Landau 0^ double-well potential. Their surprising result was the uncovering 
of a type of second-order phase transition in activation behavior at a critical value Lc of 
the system size. That a crossover in activation behavior must take place is clear from both 
simple physical and mathematical arguments (cf. Sec. 2). What is not so obvious is that 
the crossover is an asymptotically sharp, second-order phase transition in the limit of low 
noise. The change of behavior arises from a bifurcation of the transition state, from a zero- 
dimensional (i.e., constant) configuration below Lc, to a spatially varying (degenerate) pair 
of "periodic instantons" above Lc- 

The quantitative effects of the transition are significant. In the weak-noise limit, the 
activation rate is given by the Kramers formula F ~ Fq exp(— Al^/e), where e is the noise 
strength, Al^ the activation barrier, and Fq the rate prefactor. The barrier AW is inter- 
preted as the height, in dimensionless energy units, of the transition state, and by analogy 
with chemical kinetics, the exponential falloff of the rate is often called 'Arrhenius behavior'. 
The dependence on system size L of AW changes quahtatively at Lc- Also, the rate prefactor 
Fq diverges as Lc is approached both from above and below. Precisely at Lc, Fq becomes 
e-dependent in such a way that it diverges as e — > 0. This is 'non-Arrhenius' behavior. (For 
boundary conditions that give rise to a zero mode, such as periodic, there is in addition a 
noise dependence of Fq above Lc, and the prefactor divergence as L — > may be affected.) 

Given the increasingly anomalous behavior of the escape rate as Lc is approached from 
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either side, a few words should be said about the domain of vahdity of the Kramers formula 
r ~ roexp(— Al^/e) for the escape rate, which displays Arrhenius behavior with both Fq 
and AW independent of e. Strictly speaking, this formula is asymptotically valid, that 
is, only in the limit e — > 0. In a looser sense, the formula can often be applied to physical 
situations when the noise strength e is both small compared to AW, and so that the pref actor 
Fq is small compared to exp{—AW/e). These represent minimal requirements; in all cases 
applications need to be made with care. For a fuller discussion on these and related issues, 
see [21]. For the models presented here, a discussion of the regions of vahdity of all derived 
formulae will be presented in Sec. 4.3. 

A question that naturally arises is whether this critical behavior is generic: could it de- 
pend on special features of the potential studied in [20], in particular its — > —(f) symmetry? 
It was noted [^^1 that in more complicated models'-*^^] , the transition may become first-order; 
in others, it could conceivably disappear altogether. The purpose of this paper, however, is 
to provide support for the claim that the transition found in [20] is at least not confined to 
models with — > —(p symmetry; that it should in fact appear in a wide range of models 
and corresponding physical situations. To support this, a nonsymmetric (p^ model will be 
studied and solved, and a second-order transition similar to that found in [20] will be uncov- 
ered. This will be followed by a brief discussion of general nonsymmetric Ginzburg-Landau 
models with smooth polynomial potentials up to degree four, and it will be argued that this 
second-order transition should appear in typical representations of these models. 



2 The Model 

We consider on [— L/2,L/2] a classical field t) subject to the potential 

= + ^70' (1) 

as shown in Fig. 1. 

The time evolution of the field is governed by the stochastic, overdamped Ginzburg- 
Landau equation 

- Kcp" + a--f(l)^ + e^/^^{x, t), (2) 

where ^(x,t) is white noise, satisfying (^(xi, ti)^(a;2, ^2)) = S{xi—X2)S{ti—t2). The zero-noise 
dynamics satisfy = —6T-C/6(j), with the energy functional 



r-L/2 

L/2 



dz. (3) 
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V(u) 




Figure 1: Potential for the reduced field u{x,t). 



Scaling out the various constants by introducing the variables x — (0:7) ^/^/«;^/^ 



z, u 



yf^(t>, and Eq = K^l'^ot'l^l-f'l^ yields 
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2'" 



/\2 



dx 



(4) 



where £ 



L. 



It is already clear that a crossover in activation behavior must occur. In the limit £ ^ 
the gradient term in the integrand of the energy in Eq. (4) will diverge for a nonuniform state; 
while for £ ^ 00 the V{u) term will diverge for a uniform state. In this paper we will employ 
periodic boundary conditions throughoTit. and it is clear that there must be a crossover from 
a uniform to a nonuniform transition state as £ increases from 0. Physically, the crossover 
arises from a competition between the bending and bulk energies of the transition state. 

This crossover will be analyzed in succeeding sections; we will see that it corresponds to an 
asymptotically sharp, second-order phase transition in the activation rate. Both stable and 
transition states are time-independent solutions of the zero-noise Ginzburg-Landau equation, 
that is, they are extremal states of 7i[(f)], satisfying 



u" = -1 



u 



(5) 
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As already noted, we will assume periodic boundary conditions throughout. So there is 
a uniform stable state Us — +1, and a uniform unstable state Uu — —1- In the next section 
we will see that the latter is the transition state for i < i^ — v^tt. At 4 a transition occurs, 
and above it the transition state is nonuniform. 



3 The Transition State 

Following the notation of [20], we denote by Uinst,m{x) the spatially varying, time- independent 
solution ("instanton state") to the zero- noise extremum condition Eq. (5), for any m in the 
range < m < 1. The instanton state is (see Fig. 2) 



/ \ _ (2 - 3 2 

Uinst,m[X,Xo) — == / g f^^^ 

V — m+ 1 ym'^ — m-\-l 



X — Xq 



V2(m2-m + l)V4 



m 



(6) 



where dn(-|m) is the Jacobi elliptic function with parameter m, whose half-period equals 
K(m), the complete eUiptic integral of the first kind^^^l. Accordingly, imposition of the 
periodic boundary condition yields a relation between I and m: 

£ = 2V2(m2-m + l)^/^K(m). (7) 

The minimum length that can accommodate this condition is £c = a/27i", corresponding 
to m = 0. In this limit, dn(a;|0) = 1, and the instanton state reduces to the uniform unstable 
state M„ = — 1. As m — > 1~, £ — > cxD, and the instanton state becomes 

^ j (8) 

as shown in Fig. 3. As we will see, the instanton given by Eq. (6) is a saddle, or transition 
state, above Ic] it will be seen to have one unstable direction (in addition to a zero mode 
resulting from translational symmetry) . Physically, it can be thought of as a pair of domain 
walls, each of which separates the two uniform states over a region of finite extent. 

For small but nonzero noise strength the leading-order asymptotics in the escape rate from 
the metastable well are governed by the energy difference between the transition state, which 
is — — 1 below (,(. and liinst ,m{x) above, and the stable state Ug — ^-l- (A stability analysis 
justifying these identifications will be given in Sees. 4.1 and 4.2.) In the Kramers formula in 
Sec. 1, the activation barrier Al^, which governs exponential dependence of the escape rate 
on the noise strength, equals twice this energy difference. So below 4, ^W/2Eq — (4/3)^. 
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Figure 4: The activation barrier AW/2Eq for periodic boundary conditions (solid line). 
The dashed line indicates the crossover from the uniform transition state to the instanton 
transition state at £c, and the dotted line an extension to the region i > of the activation 
energy corresponding to the uniform state. 



Above £c, we find 
AW 



2En 



= ('H[Uinst,m{x)] - H[Us = 1]) 



12V2 



5(m2-m + 1)1/4 



2 - 3m - Sm^ + 2m^ 2 
3(m2 — m + 1)^/^ 3 

, (2 — m)(l — m) , , 
2E(m) - \ „ TT^K(m) 



(9) 
(10) 



(m^ — m + 1) 

where E(m) is the complete elliptic integral of the second kind^^^l. The activation barrier 
for the entire range of £ is shown in Fig. 4. As £ — > oo, AW/2Eo — > 24-\/2/5; this value is 
simply the energy of a pair of domain walls. 



4 Rate prefactor 

For an overdamped multidimensional system driven by white noise, the rate prefactor Fq 
can be computed as followsl^-*^' ^"^1 (see also [2, 14]). As in [20], let (ps denote the stable state, 
and let denote the transition state; it will be assumed (as is the case here) that this state 
has a single unstable direction. Consider a small perturbation 77 about the stable state, i.e., 
(f — (fis + r). Then to leading order 77 = — A^Ty, where A^ is the linearized zero- noise dynamics 
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at ips- Similarly is the linearized zero-noise dynamics around Thent^-*^' ^"^1 



1 



27r\ 



det A, 



I A. 



ti,i I ) 



(11) 



det A^ 

where A„^i is the only negative eigenvalue of A„, corresponding to the direction along which 
the optimal escape trajectory approaches the transition state. In general, the determinants 
in the numerator and denominator of Eq. (11) separately diverge: they are products of an 
infinite number of eigenvalues with magnitude greater than one. However, their ratio, which 
can be interpreted as the limit of a product of individual eigenvalue quotients, is finite. 



4.1 f < 4 

In this regime, both the stable and transition states are uniform, allowing for a straight- 
forward determination of Fq by direct computation of the determinants. Using reduced 
variables, the stable state Ug — -1-1 and the transition state Uu — — 1- Linearizing around the 
stable state gives 

f] = -AM77 = -{-dF/dx" + 2)77 , (12) 

and similarly 

f]^ -K\uu]r] = -{-d^/dx^ -2)-q (13) 
about the transition state. The spectrum of eigenvalues corresponding to K\u^ is 



A! = 2 + 



£2 



and the eigenvalues corresponding to h.[uy\ are 



K = -2 + 



^2 



n 



n 



0,±1,±2. 



0,±1,±2.. 



(14) 



(15) 



This simple linear stability analysis justifies the claims that Ug is a stable state and Uu a 
transition state, or saddle point. Over the interval [0, £c) ^ eigenvalues of K[u^ are positive, 
while all but one of h[uy\ are. Its single negative eigenvalue Ag = —2 is independent of I, and 
the corresponding eigenf unction, which is spatially uniform, is the direction in configuration 
space along which the optimal escape path approaches u^. 

Putting everything together, we find 



To = - 

TT 



1 If 



(2 + 



1 Ir: 



1 sinh(£/^) 
TT sin(£/V2) ' 



-2 + 



P 



(16) 



8 



which diverges at Ic — V^n as expected; in this hmit, Fq ~ const x {ic — l)"^- The diver- 
gence arises from the vanishing of the pair of eigenvalues A^^ as £ — > £j (each eigenvalue 
contributing a factor {ic — This indicates the appearance of a pair of soft modes, 

resulting in a transversal instability of the optimal escape trajectory as the saddle point is 
approached. 



4.2 i>ic 

Computation of the determinant quotient in Eq. (11) is less straightforward when the tran- 
sition state is nonconstant. This occurs when i > ic, where the transition state Uu is given 
by Eq. (6), and its associated linearized evolution operator is 



A[„ ] - + 2(2 - m) 6 



X — Xq 



A/2(m2-m + l)V4 



m 



(17) 



Evaluation of Fq therefore requires determination of the eigenvalue spectrum of A[mu] with 
periodic boundary conditions. 

An additional complication follows from the infinite translational degeneracy of the in- 
stanton state (i.e., invariance with respect to choice of xq). This implies a soft collective mode 
in the linearized dynamical operator A[-u„] of Eq. (17), resulting in a zero eigenvalue. Removal 
of this zero eigenvalue can be achieved with the McKane-Tarlie regularization procedure'^^^ 
for functional determinants. 

That procedure is implemented as follows (see [18] for details). Let yi{x,XQ;m) and 
y2{x,Xo;m) denote two linearly independent solutions of A[uu]yi = 0, i = 1,2. Let det'A 
refer to the functional determinant of the operator A with the zero eigenvalue removed. 
Then, with periodic boundary conditions, it is formally the case that 

det^A ^ y^{z + i)-y2{z) 

{yM y,iz)detlliz) ^ ^' 

where z is arbitrary, {yi\yi) = /_f(/2 '^^ Vii^t 0; m) is the square of the norm of the zero mode 
and detH(2;) = y2{z)yi{z) — yi{z)y2{z) is the Wronskian. (The expression (18) is meaningful 
only as part of a determinant quotient; see below.) 

Solutions yi and y2 of K[uy\yi = can be found by differentiating the instanton solution 
(6) with respect to xq and m, respectively; i.e., yi{x,XQ;m) = duinst,m{x,XQ]m)/dxQ and 
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y2{x,XQ]m) — dui^st,mix,xo;m)/d'm. The homogeneous solutions are then 



yi{x,Xo;m) = 
and 

y2{x,Xo;m) -- 



6m/3(m)" 



-sn 



P{m){x — Xq) 



V2 



m 



cn 



P{m){x — Xq) 



V2 



m 



dn 



P{m){x — Xq) 



V2 



m 
(19) 



3m/3(m)6 3/3(m) 



+ 



(1 -m) 



sn 



(5{m){x - Xq) 



^ 3(2m-l)/3(m)^ ^^, 



(3{m){x - Xq) 



m 



f3{m){x - xo) 



V2 



m 



V2 



+ 3p{m)hn 



P{m){x — Xq) 



V2 







m 


cn 







P{m){x — Xq) 



V2 



m 



dn 



(3{m){x - Xq) 



V2 



m 



X 



2V2 



1 — m 



(20) 



where (3{m) = (m^ — m + 1)-!/^ and E( ■|m) is the incomplete elliptic integral of the second 
kind^^^l. Inserting these solutions into Eq. (18) yields 



det' Ab, 



{yi\yi) 



(m 



2 1)11/4 



2E(m) (2 - m)K(m) 



1 — m w? — m + 1 



(21) 



9m^(l — m) 

Using a similar procedure (see the Appendix of [18]), we find the corresponding numerator 
for the determinant ratio in Eq. (11) to be 



detA['u,] = 4sinh2(£/\/2) , 



(22) 



consistent with the numerator of Eq. (16), obtained through direct computation of the 
eigenvalue spectrum. We emphasize again, however (cf. the discussion above Eq. (16)), that 
it is only the ratio of the determinants that is sensible, not the individual determinants 
themselves: these each diverge for every I. In contrast, the expressions in Eqs. (21) and (22) 
are well-behaved for all finite £ > ("^ > 0), but still separately diverge in the £ — > 00 
(m — > 1) limit. Nevertheless, here also the divergences cancel to give 



lim 

m— >1 



dct Af-u, 



24v^/J(z/i 



\Vi) 



(23) 



det' K[uy\ 

We next compute the eigenvalue A^^i corresponding to the unstable direction. With the 
substitution w — (3{m){x — Xo)/y/2, the eigenvalue equation A[ii„]77 = Xr] becomes 

d^r]/dw' + 12dn2[ ■u;|m]?7 = Srj (24) 
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where E — —2\/(5{mY + 4(2 — m). Using the identity^^^l dn^[2;|m] — 1 — msv?[z\m\, we 
observe that Eq. (24) is the I — 3 Lame equation t^^' , a Schrodinger equation with periodic 
potential of period 2K(m). Its Bloch wave spectrum consists of four energy bands, and 
its eigenfunctions can be expressed in terms of Lame polynomials^'^^\ A fuller discussion, 
especially for higher Z-values, is given in [26] ; for our purposes here we do not need to utilize 
the full machinery of the Hermite solution (a detailed treatment is given in [24]). 
It is easy to check that the eigenvector r^^^ with smallest eigenvalue is 



■nu,\{x,XQ\m) = dn^ 



(3{m){x — Xq] 



V2 



m 



+ C(m) dn 



P{m){x — Xq) 



V2 



m 



(25) 



where C(m) = -[2(2 - m)/5] + (l/5)V4m2 - m + 1, and 

1 



A, 



(2 — m) + V Aw? — m + 1 



(3{mf 



(26) 



which approaches —2 as m ^ in agreement with the single negative eigenvalue of Eq. (15). 

As noted in [20], the full translational symmetry of the periodic instanton state in the 
periodic boundary condition case corresponds to a 'soft mode', resulting in appearance of 
a zero eigenvalue \u,2 = of the operator A[itu]. (The corresponding eigenfunction 77^ 2 
is given by yi in Eq. (19).) Physically, this corresponds to zero translational energy of 
the instanton; i.e., it can appear anywhere in the volume (in contrast with, say, Dirichlet 
boundary conditions, where the instanton is 'pinned'.) This should result in an overall factor 
of £, indicating that the physical quantity of interest above £c is the transition rate per unit 
length. The general procedure for including this correction is described by Schulman [14]. 
(Our case differs by a factor of 2 from his due to the lack of symmetry in our model.) The 
net result is to multiply the prefactor by 



(27) 



The most important qualitative changes are the e~^/^ factor, leading to a non-Arrhenius 
transition rate above ic, and the effect on the behavior as £ — both will be discussed in 
more detail below. 

The above discussion also makes clear that it is not necessary to separately evaluate 
{yi\yi). For completeness' sake, however, we present it as well. Its evaluation is straightfor- 
ward and yields 



36V2(3{mf{m^ - 3m + 2) 



15 



15 



K(m) 



(28) 
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Figure 5: The Kramers rate prefactor Tq. Because of the quahtative change in behavior, 
from Arrhenius to non-Arrhenius, at ic, different vertical scales are used below and above 
the transition. 

which equals 24-\/2/5 in the m — > 1 hmit. 

Putting everything together, we find the rate per unit length prefactor for i > £cto be 



The prefactor over the entire range of £ is plotted in Fig. 5. The prefactor divergence has 
a critical exponent of 1 as £ ^ i~. Above £c the prefactor is non-Arrhenius everywhere, and 
the vertical axis is rescaled to account for the singular e^/^ behavior. 

The rescaled prefactor above £c as a function of m is shown in Fig. 6, in order to indicate 
more clearly the m — > — > £c) behavior. 

The behavior of the rate prefactor Fq for all £ > £c is unusual in two ways. First, it is 
non-Arrhenius — that is, it scales as e"^^^ for all e ^ 0. Second, it does not formally diverge 
as £ — > as seen in Fig. 6 (in fact, the divergence is present but 'masked', as discussed 
below) . Both of these features are consequences of the translation- invariance of the periodic 
boundary conditions used here, and would not appear if translation-noninvariant boundary 
conditions, such as Dirichlet or Neumann, are used. (See, for example. Fig. 3 of [20].) 



ro/£ 



X 




(29) 



\ 2(m2 -m + l)E(m) - (2 - m)(l - m)K(m) ' 
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Figure 6: The Kramers rate pref actor Tq per unit length above £c- 



Such boundary condition-dependent behavior should be distinguished from the boundary 
condition-independent formal divergence of the prefactor as the critical length is approached, 
as explained below, and of the non-Arrhenius behavior (to be discussed in Sec. 4.3 below) 
exactly at the critical length. By boundary condition- independent, we mean behavior that 
is seen in all four of the most commonly used boundary conditions in this type of problem, 
namely periodic, antiperiodic, Dirichlet, and Neumann; all were considered for symmetric 
quartic potentials in [27]. In the present case of periodic boundary conditions, the removal 
of the zero mode that is present for all i > ic renormalizes the prefactor by the factor ( in 
Eq. (27). This renormalization masks the divergence of the determinant ratio as £ — >■ 
because the factor ( includes the Jacobian of the transformation^^^] from the translation- 



invariant normal mode to the variable x^; this in turn equals the norm \J{yi\yi), which 
vanishes as m ^ 0. The crucial point is that a divergence is still embedded within the pref- 
actor, in the sense that the square root of the determinant ratio diverges with a critical 
exponent of 1/2 as £ ^ Upon closer examination, this arises from the lowest stable 
eigenvalue, A„ 3, of det A'[m„] approaching zero as f in a similar fashion to the eigen- 

value behavior below £c (cf- n = ±1 in Eq. (15)). This eigenvalue and its corresponding 
eigenfunction r]u,z are given by 




'nu,2.{.x,XQ]m) = dn' 



3 \ (3{m){x - xq) 



m -\- C"(m)dn 



(3{m){x - Xq) 
V2 



m 



(30) 
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Figure 7: Behavior of the lowest three eigenvalues of the operators A[?i,J^] and A[-u^]. The 
vertical dashed line is £ = Ic- The eigenvalues correspond to n = 0, ±1 of Eq. (15). 
Computation of the three lowest eigenvalues \i2z in the i > £c region is given in the text. 
The bold curve corresponding to A^^ reflects its double degeneracy. 



where C"(m) = -[2(2 - m)/5 + (l/5)V4m2 - m + 1], and 

"1 



A, 



u,3 



(2 — m) — V 4m^ — m + 1 



(3{mf . (31) 



Fig. 7 shows the lowest three eigenvalues for the operators A[m^] and A[m^], where 
{u^) indicates the transition state below (above) l^- This figure illustrates the evolution of 
the eigenvalues (A^j^ and A3 ) that control the formal prefactor divergence as £ passes through 
£c- We note in particular the merging of the second and third eigenvalues of A[m^] as £ — > 
consistent with the double degeneracy of the corresponding eigenfunction when I < l^.. The 
eigenvalues are everywhere continuous. Fig. 8 displays the behavior of the full determinant 
ratio above £c- 



4.3 Interpretation of the Prefactor Divergence 

The formal divergence of the Kramers rate prefactor at a critical length Ic (cf. Fig. 5) re- 
quires interpretation. It is interesting that a prefactor divergence was also found^^^' in 
a completely different set of systems, namely spatially homogeneous (i.e., zero-dimensional) 
systems out of equilibrium, i.e., in which detailed balance is not satisfied in the stationary 
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Figure 8: Behavior of the ratio of the eigenvalue spectrum above Re- 
state. That divergence arose from an entirely different reason: the appearance of a caustic 
singularity in the vicinity of the most probable exit path as a parameter in the drift field 
varied. The caustic singularity arises from the unfolding of a boundary catastrophe; a de- 
tailed analysis in given in [31]. In contrast, the problem considered here is that of a spatially 
extended system in equilibrium, so no such singularities can be present. Moreover, no pa- 
rameter in the stochastic differential equations describing the time evolution of the system 
is being varied in the case under discussion here; rather, the variation is in the length of the 
interval on which the field is defined. The 'phase transitions' in the stochastic exit problem 
in the two classes of systems are therefore physically unrelated. 

What does it mean for the prefactor to (formally) diverge? In fact, at no lengthscale is 
the true prefactor infinite, for any e > 0. Indeed, given that the analysis presented here is, 
strictly speaking, valid only asymptotically as e — > 0, the escape rate is always small where 
the above results are applicable. What the formal divergence of the prefactor does mean 
is that the escape behavior becomes increasingly anomalous as £c is approached, and it is 
asymptotically non-Arrhenius exactly at 4- 

That is, when £ = ic the true rate prefactor should scale as a (negative) power of e 
for all e ^ 0. As in [28, 29], this can be treated quantitatively by studying the 'splayout' 
of the 'tube' within which fiuctuations arc largely confined, as the saddle is approached. 
'Splayout' here simply means that the fluctuational tube width, which for £ 7^ £c is 0(e^/^), 
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becomes 0(e"), with a < 1/2, as the saddle is approached when £ — £c. In the model 
studied in [28, 29], the Lagrangian manifold comprising optimal fluctuational trajectories 
has a more complicated behavior than in the model under study here. As a result, in [28, 29] 
fluctuations near the saddle occur on all lengthscales, while in the present case fluctuations 
near the saddle occur on a definite lengthscale, but one larger than 0(e^/^), leading to non- 
Arrhenius behavior for all e — > 0. Details will be presented in [30]. 

Our main interest here is in the region close to ic, where the rate pref actor Fq is growing 
anomalously large (but remains everywhere finite for all i strictly away from ic). The for- 
mulas given here are valid for e sufficiently small so that the contribution from the quadratic 
fluctuations about the relevant extremal state of 7Y[0] dominates the action. As long as all 
eigenvalues of A[uy\ are nonzero (excluding the zero mode arising from translational sym- 
metry, which may be extracted), and the norms of the corresponding eigenfunctions are 
bounded away from zero, the pref actor formula applies, but in an e-region driven to zero as 
£ — > by the rate of vanishing of the eigenvalue(s) of smallest magnitude. Therefore, as 
£ ^ ic from either side, the Kramers rate formula applies when e scales to zero at least as 
fast as \i — ic\^^'^ (of course, the constraints already mentioned in Sec. 1 must continue to 
hold as well) . More precisely, the criterion considered here (which is necessary but not a pri- 
ori sufficient) for Arrhenius behavior to hold on either side of £c is that the noise strength 
e be small compared to Xm{Vm\Vm) , where is the eigenvalue of smallest magnitude and 
Tjm its corresponding eigenfunction(s). This criterion arises from the condition, used in the 
derivation of Eq. (11), that the noise strength be small compared to the size of quadratic 
fluctuations about the extremal action. For i slightly below ic, these quantities are given in 
Sec. 4.1, and above 4, by Eqs. (30) and (31). The resulting computation is straightforward 
and the resulting e-region is sketched in Fig. 9 (in the dimensionless units used here, the 
coefficients of the scaling terms are of 0{1)). 

The result is that, for fixed i close to ic, there should be a crossover from non- Arrhenius 
to Arrhenius behavior at sufficiently weak noise strength (cf. [29]). The figure represents a 
type of 'Ginzburg criterion' that describes, at a given i near ic, how far down as e ^ the 
non- Arrhenius behavior persists. It should be emphasized that Fig. 9 sets an upper bound 
on the scaling of the region of e vs. \i — id, below which asymptotic Arrhenius behavior sets 
in. It would be interesting to consider also the behavior at very small but fixed noise as i 
increases thorugh ic- Here one would observe a crossover from Arrhenius to non- Arrhenius 
behavior and back again as i passes through the critical region. An interesting problem for 
future consideration is to analyze this phenomenon in greater quantitative detail. 
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Figure 9: A sketch of the scahng of the regions where the Arrhenius prefactor formulae given 
by Eqs. (16) and (29) are vahd, when I is very close to £c- For fixed £, e must be small enough 
so that it is below the shaded region for these formulae to apply. Non- Arrhenius behavior 
is expected in the shaded region. Because this figure is intended to illustrate only the rate 
of scaling of e with | / — so that the Kramers rate formula is valid, the axes are unmarked 
(except for £ = where the e-range has shrunk to zero). 

To summarize: strictly away from the prefactor formulae Eqs. (16) and (29) hold (cor- 
responding to Arrhenius behavior of the rate), but in an increasingly narrow range of e as 
£c is approached. A crossover from non- Arrhenius to Arrhenius behavior as e — >^ should be 
observed, along a boundary that scales as shown in Fig. 9. Strictly at i^^ the formulae do not 
hold: the prefactor is finite, but acquires a power-law (in e) character. That is, the rate be- 
havior is non- Arrhenius all the way down to e ^ 0. This (boundary condition- independent) 
non-Arrhenius behavior at criticality should be distinguished from the (noncritical) non- 
Arrhenius behavior strictly above £c that appears only when translation-invariant boundary 
conditions, such as periodic, are used. 

5 Asymmetric quartic potentials 

In this section we will consider more general asymmetric potentials of the form 



V{u) 



(l/2)«2 _ (i/3)«3 + (1/4)^4 



(32) 
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Figure 10: The asymmetric quartic potential of Eq. (32). 



as shown in Fig. 10. 

We will consider only the small-£ regime, and show that a transition at finite £c exists 
with a divergence of the prefactor as £ — > £c- 

Because the prefactor depends on the curvature of the potential near its minimum, dif- 
ferent prefactors (and of course activation energies) correspond to the two minima u_ = 
1/2 — \/5/2 ^ —.62 and u+ = 1/2 + \/5/2 ^ 1.62. The formulae shown here correspond to 
U-] corresponding formulae for escape from the potential minimum from u+ are obtained by 
replacing the constant a- — \1 + 2m_ — 3m^| fti 1.38 with a+ — \1 + 2u+ — 3m^| fti 3.62. The 
critical length £c is the same in the two cases. 

An analysis similar to that of Sec. 4.1 yields 



so £c — 27r and, as before, the prefactor again diverges as {1 — i/ic) ^ as £ — > , as shown 
in Fig. 11. 

6 Conclusion 

We have found an explicit solution of the Kramers escape rate in an asymmetric (f)^' field 
theory of the Ginzburg-Landau form. This result, and the brief discussion in Sec. 5 of 



ro 



(a_)V^ sinh(a_£/2) 
sin(£/2) 



(33) 
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Figure 11: The prefactor Tq of Eq. (33) for £ < 4. 



more general asymmetric potentials, suggests that the critical behavior found in [20] might 
hold for a more general class of models than those with a high degree of symmetry. How 
widespread the transition phenomenon is remains uncertain, but it appears to hold at least 
for arbitrary smooth potentials with terms up to and inchiding 0^. It would be interesting 
to find models with other types of behaviors. One interesting possibility, discussed in [20], 
is a class of models that display a first-oideT transition: for example, a discontinuity in the 
derivative of the activation barrier height with respect to the interval length, at a critical 
length. A possible candidate for such a model is the sixth-degree Ginzburg-Landau potential 
of Kuznetsov and Tinyakov [19], but a detailed analysis of its transition behavior remains 
to be done. 
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